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ABSTRACT 

Extensions of the investigations aimed at preventing aircraft destruc- 
tion through the effects of hydraulic ram were conducted. By analyzing 
spark shadowgraphs the existence of a family of transient shock waves 
was confirmed. These shock waves were created by a surface anomaly in 
a projectile and are not unlike those found in fragmentary warheads. In an 
attempt to measure wall deflections subsequent to projectile impact, a 
system based on optical fiber proximity probes was investigated. Mathe- 
matical models based on grid distortions were used as a means of calculat- 
ing shock wave pressures. The AFIT thesis correlating cavity size, energy 
and radius to underwater explosions was verified from an analysis of 16 mm 
films . 
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I. INTRODUCTION 



The term hydraulic ram has been historically associated with the 
various waves causing, usually, catastrophic failure of an aircraft fuel 
cell when penetrated by a projectile. 

When a projectile or fragment with high velocity impacts the surface 
of a liquid medium, a disturbance is produced which spreads throughout 
the medium. The disturbance is characterized by a region of increased 
pressure and density together with a small forward motion of the fluid. 
When the energy of the disturbance is large, the disturbance travels as a 
surface of discontinuity and is called a shock wave. (Ref. 1). Typically 
the shock wave phenomenon has been identified as the first phase of hy- 
draulic ram. 

As projectile penetration continues, a pressure field is created that 
tends to displace the fluid from the projectile, transforming projectile 
energy to kinetic energy of the fuel. This is the beginning of the cavity 
phase. 

The phenomenon associated with these two phases has been studied 
for as long as three wars i.e. 1947 through present. Recent investigations 
into fragment-induced shock waves and cavities have been initiated, the 
results of which are leading to a more clear understanding of hydraulic 
ram. 
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II. SHOCK PHASE EXPERIMENTS 



A. WAVE GENERATION IN SHOCK PHASE BY FRAGMENTARY PROJECTILES 

Previous experiments conducted at NPS involving the impacting of 

0.22 caliber projectiles into a water-filled tank produced shadowgraphs 
such as those shown as Figures 1 and 2. Figure 1 is the projectile in free 
flight , whereas Figure 2 shows the family of waves created after impact. 
The flow features visible in Figure 2 are considered to be: 

1. projectile wake, 

2. reflected bandwave, 

3. reflected bow wave, 

4. incident wake recompression shock, 

5. incident bandwave, 

6. transient entry shock, 

7. secondary transient shock wave (see discussion below) . 

In Figure 2, items 6 to 9 occur in the water internal to the tank and items 
1 to 5 occur in the air external to the tank. The dense, thick region in 
each shadowgraph is the test tank wall. Region 8 is the projectile's 
cavity and item 9 is the grid distortion used to determine the medium's 
density. The spark shadowgraph method was used as it has the advantage 
of showing the projectile-induced waves at one observation. Of particular 
interest was the formation of the wave 7, in Figure 2, which was investi- 
gated as follows . 
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Fig. 1 0.22 Caliber Projectile in Free Flight 

1. Bow Shock Wave 

2. Band Shock Wave 

3. Wake Recompression Shock 
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Fig. 2 Waves Generated From Projectile Impact 
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A randomly selected, factory loaded 0.22 cal projectile has the geom- 
etry shown in Figure 3. The volume of the projectile's nose was calculated 
using successive truncated cones while the portion aft of the shoulder was 
calculated assuming a solid cylinder. Assuming a constant velocity immedi 
ately after penetration into the tank, the rate of volume penetration was 
plotted against distance along the projectile’s centerline , Figure 4. The 
discontinuity in Figure 4 appears at a longitudinal distance of 0.20 inches 
and corresponds to the distance where the projectile itself has a discontin- 
uity, i.e. , the shoulder. This shoulder is the origin of bandwave shown 
in Figure 1. This suggests that although the surface irregularity is small 
(0.020 inches), it is still capable of creating secondary transient shock 
waves. It is concluded that the wave 7 in question on Figure 2 is indeed 
a secondary transient entry shock. Reference 2 describes the deformation 
of several previous test rounds. Considering their fragmentary shape it 
could be expected that fuel cells punctured by fragmentation warheads 
would be subjected to a myriad of such secondary shocks . 

B. PRESSURE BY GRID DISTORTION 

Figure 5 is a shadowgraph of a fine-mesh grid taken while the fluid 
was at rest. No grid distortions are apparent. The projectile's path is 
from the left to right through the center of the shadowgraph. After penetra- 
tion the grid in this area is distorted. This region was used in the subse- 
quent analysis . 

As the projectile penetrates the water it displaces fluid by compressing 
it to very high pressures. This disturbance is the transient shock 6 in 
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Time (Microsec) 

Fig. 4 Volume Penetration vs Time for Projectile 




Fig. 5 Fine Mesh Grid Showing Planar Regions 
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Figure 2. As penetration continues, the transient shock continues to receive 



energy until the projectile becomes fully wetted; this occurs at 23 micro- 
seconds in Figure 4. A rarefaction wave then forms that tends to dissipate 
the shock wave. (Ref. 3). The localized high shock pressures are known 
to be of short duration, usually tens of microseconds, and are completely 
dissipated by the time significant projectile penetration occurs. 

The spark shadow-graph could be set off with any time delay up to 
approximately 30 microseconds. Each shock wave in the shadowgraph re- 
veals a distinct presence-of-light band and an absence-of-light band. The 
concentration of light around the projectile's transient entry shock is best 
explained by considering the pressure and density increase due to the shock 
wave. Similar optical distortion at the nose of the fragment is due to high 
density at the nose. See Figure 2. The pressure is known to reach values 

r\ 

in excess of 4000 kg/cm z (56,890 psi) , which is sufficient to produce a 
change in the medium's refractive index. The width of the dark band varies 
with the pressure in the wave. The fractional change in refractive index, 

, has been studied by Poindexter and Rosen and was determined to vary 
as € = 1.06 x 10"°p where p is pressure in kg/cm (Ref. 1) . This conforms 
to an isothermal change, whereas a change produced by the transient shock 
is more nearly adiabatic (Ref. 1) . A value of £,= 10”^p has been justified by 
McMillen and was used throughout this analysis (Ref. 1). It was also as- 
sumed that in a homogeneous medium such as water, light rays have the 
form of straight lines. 
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Fig. 6 Shock Front Pressure Generated in Water by Projectile Impact. 
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Seven shock pressures were extracted from Figure 6. When the pro- 
jectile had penetrated three cm into the tank the pressure was noted as 
2700 kg/cm^. Define©, and Q^as the angles which the incident ray and 
refracted ray make with a normal to the shock wave surface (Fig. 7). 




Fig. 7 Light rays in homogeneous regions 
bounded by a shock wave. 



By geometry is given by (1) 

sm — 

e «= sm "(V) 

Using Snell's law the angle ©-j^can be introduced 

I 

n^sw.9^- 
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From the previous paragraph the change in refractive index can be expressed 
as 6 s 

( 2 ) 
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From equations (1) and (2) 
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Note that the changes in angle entering and leaving the shock wave are 
additive. If the density is homogeneous within the region bounded by the 

shock wave, as is assumed here, then the high pressure region behaves 

P 

as a spherical lens. For-j^ =0.90 from equation (1) , ■6 J is found to be 
5.740°. For a pressure of 2700 kg/cm^, -©^is calculated to be 5.588° by 
equation (3); consequently 1LQ = 0.30 3°. Even for high pressures the 
angular deflections are small. Calculated values for the remaining six 
pressures are presented in Table 1. These tabulated results appear reason- 
able . 

The grid distortion analysis assumes rays have spherical symmetry 
i.e. the refractive index depends only on the distance r from a fixed point 
0; see Figure 8. 




For rays situated in a plane through the origin, nrsin^=c, where c is 
a constant. With d = rsin^ , nd = c 

This is known as the formula of Bouguer. From geometry, if (r,-0 ) are 
the polar coordinates of a plane curve, the angle (p between the radius 
vector to a point P on the same curve and a tangent at P is given by 
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Table 1 


Angular deflection of light ray by a spherical shock wave 
for various shock wave pressures 
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Table 1 Angular deflection of light ray by a spherical shock wave 
for various shock wave pressures 



RR Nondimensional radius 
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Table 1 Angular deflection of light ray by a spherical shock wave 
for various shock wave pressures 

RR Nondimensional radius 

Thetal Entry angle, radians 

Theta2 Exit angle, radians 

Ansi Entry angle, degrees 

Ans2 Exit angle, degrees 
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It is possible to eliminate <j> by introducing the constant c, as can be seen 
by the following manipulations: 




/. = s' + JL jk J1 

\J&~ cF * 
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The assumption of constant pressure behind the shock wave is not 



realistic. At this point, the assumption of uniform conditions behind the 
shock wave can be abandoned by introduction of equation (5) . In equation 
(5), c is determined by the point where the ray intersects the shock wave. 

The index of refraction, n, which is related to pressure by £.= 10~^p, can 
be made a function of the radius. 

To adapt equation (5) to the specific case of the hydraulic ram shock 
phase it is necessary to carefully define the various geometric parameters 
involved. Figure 9 illustrates the angles and radii necessary to apply 
equation (5) . Note that the diagram is symmetric about r* at angle . 

A functional form for n(r) can be assumed. For example, a parabolic 
relation would be 

n(r) = ar^ + br + c 

In order to evaluate the constants a, b, and c, the integral of equation (5) 
must be evaluated and iterated with geometric relations . This is quite 
complicated . In lieu of a parabolic relation a linear formula is used. This ’ 
has the merit of a straightforward solution. One of the two parameters, i.e. 
A, can be determined leaving B to be found as the following manipulations 
show. 

Assume (Fig. 9) 

TV- r* y ' , k » T\i 

V 



1 \ 



23 



* • 



( 6 ) 



*•- 'H 4 KV 0 

From the law of Bouguer 
S UA — constant 

A^V£va^<£)^s= A? v* S tu. * = K*Y* 

(7) 

From analytic geometry and Snell's law 

= *£ (8) 

iUA^ v Y-, 



$r $£ + *$- je - g*-^-^** 

«*= UE^+'LLfy 

$,*31 -Q"-2,-a^ +1^4 1^ 

*1^ 

/. TV - 0" H 2j ■+ a£> 
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( 10 ) 
(ID 
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(13) 



From Figure 9 with *,<{0 /,* ,?H known 

•io 4 Hi 4v h-~ ' r ' s ' M ^i 

VHr r ’ s ^^ = r'lui ((?>,- *) 

l^-j = x' SWa^' * f SVM,^*"<*) “ T> - iUA^-oC 

Cfili £,= TT-*A<£ 

<j^ = TT -t a£ - £ , 

o< - n t • 

L r \ 
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Fig. 9 Geometry of Light Ray 






and 
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With the following substitutions, the solution for *0 is 
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Using a pressure behind the shock wave of 400 kg/cm^, ri 2 was 
found to be 1.335. See Figure 6 for data used in these calculations: 
r' is 7.62 cm. As shown in the Appendix the value for B is -0.04947. 

Using the equation for the linear pressure distribution it is calculated 

. o X 

that p is 3646 kg/cm^ at y.* 0.1267. The pressure increases toward 
the projectile nose. This result is consistent with Figure 29 of Holm's 
thesis (Ref. 2). 

Results from fluid mechanics indicate that linear relations for p(r) 
and n(r) are not appropriate. Consequently the calculated pressure at 
~p* # = 0.12 67 is probably not correct. This analysis gives the correct 
trends. The difference between experimental and calculated values is 
considered to result from the assumed distribution for n. 

C. WALL RESPONSE BY FIBER OPTICS PROXIMITY PROBES 

The structural response of the entrance wall during the shock phase 
has been investigated and is presented in Refs. 4,5, and 6. Since the shock 
wave is a pressure disturbance propagating through the fluid medium, it is 
also by definition a stress wave. Types of stress waves include compression, 
tension and shear waves. The compression and tension waves involve posi- 
tive and negative pressures, respectively, where the media particle motion 
is in the direction of wave propagation. The shear wave involves particle 
motion that is normal to the wave direction. Compression and tension waves 
are also known as dilatation waves while shear waves are sometimes called 
distortion waves. Tension waves that reduce compressed areas are also 
known as rarefaction waves (Ref. 3). Loads transmitted through the tank's 

1 \ 
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structural members are initially transmitted as stress waves of both dilata- 
tion and shear form. The outward motion of the entrance wall somewhat 
relieves the pressure disturbances (stress waves) acting on the wall. To 
measure the outward displacement of the wall a system based on optical 
glass fiber bundles was constructed. The optic fiber bundles were arranged 
coaxially with the transmitting bundle on the inside. This light was then 
directed onto the aluminum entrance wall. The reflected light was collected 
and transmitted by the output fiber bundle to an RCA 9 31A photomultiplier 
tube and a Pacific Photometric Model 11 laboratory photometer. The fiber 
bundles were three feet in length and were comprised of three-mil flint 
glass fibers with an index of refraction of 1.62 and a cladding glass of in- 
dex 1.52. The corresponding numerical aperture was 0.56 with an accept- 
ance cone half angle of approximately 34 degrees. The glass bundles were 
secured in a mounting bracket and placed over the entrance wall on a radial 
line emanating from the projectile impact point. The mounting bracket al- 
‘lowed for ’positioning the bundles up to a maximum of nine inches from the Ali ,, U il 
impact point. To increase the intensity of reflected light the entrance wall 
was hand polished with #500 Wetordry Tri-M-ite paper, then with Crystalbay 
Crocus cloth. The signal output of the photomultiplier was transmitted to a 
dual beam oscilloscope fitted with a Polaroid camera. 

Initial calibration of the system was accomplished by mounting the 
fiber bundles in a platform micrometer. The light was reflected from a 
circular aluminum disk of the type used for the entrance wall. This disk 
was similarly polished. A calibration curve of volts versus the distance 

separating the disk and glass bundles was obtained, Figure 10. 
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Subsequent investigation of this system, modified to include a 



"chopper wheel", revealed an unacceptably low rise time for the photo- 
multiplier tube. A Pacific Photometric Model 62 multiplier tube (RCA model 
6199) was substituted; a Spectrum Physics model 132 HeNe laser replaced 
the light bulb. The inherent noise in the latter system obliterated the os- 
cilloscope trace, rendering this optical fiber concept useless without 
further development to improve signal-to-noise ratio. 
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Fig. 10 Calibration Curve for Optical Fiber Proximity Probe 
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III. CORRELATION OF OSCILLATION PERIOD. 

CAVITY SIZE, AND ENERGY DISSIPATION DURING CAVITY PHASE 

A. TECHNICAL APPROACH TO CORRELATION 

In displacing water from its path, the projectile imparts a radial 
velocity to the water. The momentum associated with this velocity field 
is sufficiently high to cause the flow to separate from the projectile and 
to continue to expand causing a cavity. Fluid pressures caused by the 
tank boundaries eventually halt the cavity's expansion. At this point 
most of the kinetic energy initially associated with the velocity field has 
been converted to stored energy, gravitational potential energy and de- 
formation of the tank's wall. These forces cause the cavity to start to 
collapse radially and, as the cavity shrinks, the fluid velocity increases 
resulting in the conversion of stored energy into kinetic energy and work 
in compressing gas trapped in the cavity. Since projectile impact was near 
the center of the tank the cavity would first neck down until a nearly spheri- 
cal shape resulted. The sphere would then tend to collapse to a point until 
the pressures in the cavity halted this inward motion. Air which entered 
the cavity through the entrance wall hole is trapped and compressed to a 
high pressure at the instant of minimum cavity volume. A high pressure 
shock wave then propagates radially outward from the collapse point. After 
the pressure wave generated by the collapse subsides, the compressed air 
causes the cavity to re-expand. This cycle repeats several times until all 
the stored energy is dissipated (Refs. 6 and 7). Various authors, Lundstrom, 
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Stepka and Cardea, et al. , have likened the cavity phase to underwater 



explosions. Using the method presented by Cardea and Cole, an analysis 
was made of a high speed film showing 0,30 cal and 14.5 mm rounds im- 
pacting a 60 61T6 aluminum entrance wall. The film was obtained from 
Naval Weapons Center. These values were then compared with the analyti- 
cal solutions obtained by Holm's research at NPS. Two significant varia- 
tions were noted. Reference 7, equation (7), presents the maximum cavity 



accurate. 'Shis difference accounts for an 8% variation in tabulated results 
shown as Table 2. Secondly, this paper considers the effects of projectile 
impact through a pre-punched aluminum entrance wall whereas the 0.30 
cal and 14.5 mm rounds impacted a solid tank wall. For the latter the 
shock wave is generated by the impact of the projectile on the wall vice 
the fluid medium. The generation and propagation will be identical ex- ’ 
cept for wall density and sonic velocity effects. This problem is addressed 
in Reference 3. The analysis of the 14.5 mm round must be considered as . 
approximate as the maximum cavity radius to be measured extended beyond 
the frame size of 16 mm film. 

B. DATA OBTAINED FROM MOVIE FILM 

For each sequence of the 16 mm film obtained from the Naval Weapons 
Center, China Lake, the maximum cavity radius and time for one complete 
oscillation were measured. Using underwater explosion theory (Ref. 8) 



radius 




. A constant of .544 i s considered to be more 



the energy required to expand the cavity was obtained first as a function 



of cavity radius and second as a function of the cavity period. A theoret- 
ical value for the cavity radius was also calculated. These results were 
then compared with those obtained by Cardea (Ref. 7). 



C. RESULTS 

Table 2 lists the results obtained from the above measurements. The 
conclusion equating the cavity phase of hydraulic ram to underwater ex- 
plosions appears valid. 

For a comparison of the time required for one complete cavity oscilla- 
tion, Reference 2 presents 



“L .W, fyy 

rw 



where fl^is the measured cavity radius. From Reference 8 this time is 

J. 



given as 



-t= 



where 



v'* 



2210 psf 
1.04 slugs 






Evaluating the constants yields 

2 - l0 -«(kf 

t= 

respectively. It can be seen that the only difference in characteristic 
times is in the constants. Based on results of analysis of the film, 
equation is considered more correct. 
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0.30 cal 


t 


a m 


a m 


Ec=f(Am) 


Ec=f(T) 


Shot number 


millisec's 


measured 


calc 


Ft-lbs 


Ft-lbs 






Ft 


Ft 






1 


24.8 


.485 


.455 


874.78 


874.81 


2 


22.0 


.555 


.404 


610.68 


610.70 


3 


28.1 


.632 


.516 


1272.5 


1272.6 


4 


28.8 


. 656 


.529 


1370.0 


1370.1 


5 


29.1 


.632 


.534 


1413.3 


1413.3 


14 . 5 mm 
Shot number 












1 


46.9 


1.125 


.861 


5916.5 


5916.7 


2 


48.3 


1.094 


.887 


6462.3 


6462.5 


3 


47.4 


1.007 


.870 


6107.7 


6107.7, 


4 


44.4 


1.234 


.815 


5019.9 


5020.1 


5 


45.5 


1.266 


.836 


5402.3 


5402.5 


6 


43.9 


1.234 


.806 


4852.2 


4852.4 



Table 2 Cavity Period, Maximum Cavity Radius, Cavity Energy 
from Movie Film 



0 . 50 cal 


t 


a m • 


a m 


Ec 


Shot Number 


millisec's 


observed 


calc 


Ft-lbs 


2 


24.8 


.5 


.494 


876.0 


3 


26.2 


.5 


.522 


1032.0 


4 


26.9 


.5 


.536 


1117.0 


5 


26.9 


.5 


.536 


1117.0 



Table 3 Cavity Period, Maximum Cavity Radius, Cavity Energy 
from Cardea (Ref. 7) 
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VI. CONCLUSIONS AND PROPOSED RESEARCH 



A. CONCLUSIONS 

Secondary shocks will occur in fuel cells penetrated by fragmenta- 
tion warheads . 

The mathematical model proposed for analyzing the optical grid dis- 

* 

tortion yields correct trends. Further development of an assumed distri- 
bution for the index of refraction is warranted. 

Comparison of the cavity phase of hydraulic ram to underwater ex- 
plosions provides valid relationships between cavity radius , energy and 
period of oscillation. 

B. PROPOSED RESEARCH 

HIGH SPEED FRAMING CAMERA 

A high speed framing camera capable of 200,000 frames per second 
has been obtained. When incorporated into the existing test set-up, a 
complete time history of the shock front and projectile penetration will 
be available . 

STRAIN GAGES 

The existing entrance wall has been redesigned and fitted with several 
strain gages. Proposed research with this entrance wall should yield in- 
formation on wall deflections and wall bending moments caused by pro- 
jectile penetration. 



! \ 



35 



APPENDIX A 



SAMPLE CALCULATIONS FOR THE GRID DISTORTION ANALYSIS 




Fig. 12 Enlarged Grid Distortion 
Steps for obtaining a solution: 

1. evaluate equation (17) to find 

2. evaluate equations (8) and (9) to find ^ and A ^ 

3. evaluate equation (12) to get . 

4. evaluate equations (10) and (13) for © and © >f 

5. evaluate equation (14) to get 

6. evaluate equation (15) to get 

7. evaluate equation (16) for 

8. evaluate equation (11) for (p £ 

9. draw a curve of B vs r* using equation (18) 

10. using number pairs (B,r*) from step 9, evaluate equation (19) 
for the right hand side. Iterate until correct (B,r*) are found. 
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From Figs. 11 and 12 
d =0.1 inch 
L 1 = 19 . 8 inches 
L, = 9 inches 
= .35 inches 
t* = 3 ” 

P = 400 kg/cm^ 

pty - (0.1) (180) = 0.29° = .005 radians 
(19.8) 

From equation (17) 

0.35 = 3 [$ - 3 CDS ($ r +0 lW.o( 

Assume ( ^— 0^) is small 

0 . 35= 3($,-o0 - (9-3)0< 

3 £, = .35 + 9«X 

or = 0.35 + 3<^ = .1318 radians 
3 

try = 1.-55° 

— 'l,*b c b C ~ 0.1^= *7.^-C 7 
then 

o.3S = hs^n.Tio-V 9-5co5 1.2^1 (+olk o. 19) = ,34-2k 



from iteration 



Using p = 400 kg/cm^ at r' = 3 inches 

ld 5 p= o.oo^ 

^,(14.004)" l.Sb ( 1 . 004 ) 

1.3333 

.'. 1,33 _ 0.$9tO 

rv^ I.33S3 

then 

,rM1( 

.-. “?.*=> i$3° 

*(> = ^ir "7. 51°- 7.3133° 

6.0511° 

from equation (12) 

= (X - A (3 
^ 

z,- .JS- .osn = 0 . 0 ^ 

v 

from equation (10) 

9'= <)&+«, + i .0^53 + .o=>n-i.s-] 

©' = 

from equation (13) 

&*= %-*-£, + = ^I.SIS 0 

from equation (14) 

Hl = (1,51- ^9)- .033 
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from equation (15) 

•b^SiAA.I.Sv Su/- (1.^- ,V))\ 

Hv 5 

from equation (16) 

|M> + .OS-l.-bl 

<f>^= nt.'Wn'’ 

from equation (11) 

K = = " 131 - • 0£> n 

<{> “ = i.eiS3° 

An approximation for B is 

VHi 



r' r' 3 

from step (9) 




$= (- t,S333>X , 1303- .Wei) 

(,ru/i)(> .ati) 

6 - — 0 . 04 -^ 47 ° 



from the assumed distribution for n 



Y. 



K + b 




1.3335 4 0. 6+9^7 

*»V = 



when 






r 



* 



« \ 



39 



then 

t =■ * W -l— Ts. i 

*• V'* | 

e= t, 3 ns- i ,33 

1,33 

6~ 0. 034>4-(o 

but 

6= lo 

.*. p- 10^(0.03^0 

2 

p=3k4k.Ta kg/crn (calculated for a radius ratio of 
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I 

_P R- NCJt-.D.I M E'sLS I QNAL RADIUS 

XN 1 - REFRACTION INDEX OF WATER, CON STANT= 1. 33 
XN2- REFRACTIVE INDEX BASED ON SMELL'S LAW 
PPE- PRESSURES EXTRACTED FROM PREVIOUS THESIS 

IF ET A Lr_E.NI BY-ANGLE *—RA C I A NS— 

TFETA2- FXTT ANGLE, RADIANS 



ANSI- ENTRY ANGLE, DEGREES 

A A S2-* EXIT ANGLE, DEGREES 

RES-^JhETAl -IHETA2-J. I YES 2 

Cl PENS ION POF (9 ) (9) , THETA L( 9 ) , THFTA2( 9) , B( S ) , E ( 9 ) , XN2 ( 9) 

D I REN SI CM XX (9) ,ANS1 (9 ) , ANS2(9 ) ,RES (9 ) 

XN1=1.23 
CC 10 J = 1 , 9 

P E AO (5,902) PRO ( J ) — 

WRITE (6,502) P»C<J) 

„W R I TE .. ( 6 , 50 1 ) f 

SLM=0. C 

E ( J )= 1 . C F - 5 Y p R E ( ) ) 

XN2( J)=XN1*(1.0 + E ( J ) ) — 

B( J )= 1.0/ ( 1 .0 + E ( J ) ) 

CC 20 1=1,9 

SIR=SEY + .1 

PR ( I ) ---SUN 

T F E T A 1 ( I > = A P S I N ( P P ( I .) ) . . __ l 

AN S 1 ( I ) =THE T A 1 ( I ) *5/. 3 
XX(I)=e(J)*FP(I) 

TFETA2 ( T ) = APS IMXX(I) ) 

ANS21 I ) = T H E T A 2 ( I )*57.3 
RES ( I ) = ( A N S 1( I )— AMS ?. ( I ) )*2. 

WRITE (6,503) RP ( I ) , THETA1 ( I ) , T HET A2 ( I ), AN$ 1(1), ANS2( 1 > ,RES ( I ) 

20 CONTINUE 

1 0. _C C N TJ N U E 

WRITE ( fc » 50A ) XN2 

501 FCPMAT ( *0 • , AX, »RR • , 8X, • THFTA1 • , AX, ' THETA2' , AX , • ANSI' ,5X, • AN S2 ’ ,7X 

1,'PES') 

5C2 FORMAT (F 8. 2) 

5C3 FORMAT (6F10 .3 ) 

5CA FORMAT (• 0’ , AX XN2 = «,E7.4) 

STOP 
c k r 
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